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Abstract. An update on astrophysical models for nucleosynthesis via rapid neutron capture, the r 
process, is given. A neutrino-induced r process in supernova helium shells may have operated up 
to metallicities of ~ 10~ 3 times the solar value. Another r-process source, possibly neutron star 
mergers, is required for higher metallicities. 
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INTRODUCTION 

Two major mechanisms are responsible for making heavy elements beyond Fe: the slow 
(s) and rapid (r) neutron capture processes [1, 2]. This contribution focuses on recent 
developments in studies of the r process. Extensive reviews on the subject can be found 
in e.g., [3, 4, 5]. See [6] for a recent comprehensive review of ^-process studies. 

There are three major areas of r-process studies: (1) properties of extremely neutron- 
rich nuclei far from stability that are involved in the process, (2) astrophysical sites and 
associated conditions that can be used to assess the feasibility of an r process, and (3) 
frequencies of occurrences and other characteristics of potential r-process sources that 
are relevant for chemical evolution of galaxies. The second area is discussed in some 
detail below while the other two areas are mentioned only briefly. 

The basic nuclear reactions involved in the r process are neutron capture, photo- 
disintegration, and /3 decay. The crucial nuclear input for the r process includes: (1) 
neutron capture rate coefficients (v<J n v(Z,A)), which are averaged over thermal distri- 
butions of neutron and nuclei, (2) photo-disintegration rates Ay >n (Z,A), which can be 
obtained from (vo n v(Z,A — 1)) using detailed balance along with the neutron separa- 
tion energy S n (Z,A) and relevant nuclear partition functions (e.g., [3]), and (3) /3-decay 
rates Xp (Z,A) . Here (Z, A) stands for a nucleus with proton number Z and mass number 
A. As the nuclei involved are extremely neutron rich, branching ratios for emission of 
up to several neutrons following a /3 decay also play an important role, especially dur- 
ing decay back to stability subsequent to freeze-out of rapid neutron capture. If the r 
process occurs in core-collapse supernovae, neutrino-induced reactions are important as 
charged-current (CC) v e capture is equivalent to /3 decay and neutrino-induced neutron 
emission is efficient following both CC and neutral-current (NC) excitation [7]. If the 
r process produces fissile nuclei, then rates for various fission channels (e.g., sponta- 
neous, neutron-induced, /3-delayed, neutrino-induced) and the corresponding mass dis- 
tributions of fission fragments should also be included. Only a small fraction of the 



nuclei involved in the r process are accessible to experiments, so the majority of the 
above nuclear input must be calculated from theory. As this input covers a wide range of 
nuclear processes, a self-consistent theoretical framework for all the input is essential. 
Fortunately, the next generation of rare isotope beam facilities such as FRIB, RIKEN, 
and FAIR will provide crucial guidance to the theoretical efforts in addition to direct 
measurements of masses and /3 -decay rates of some key nuclei for the r process. 

By definition, an r process occurs when the neutron capture rate of a typical nucleus 
involved in the process exceeds its /3 -decay rate: 



n n (va nJ (Z,A)) > Xp(Z,A), 



(1) 



where n n is the neutron (number) density at the astrophysical site. The quantity 
(vc„,y(Z,A)} depends on the temperature T at the site while A«(Z,A) is essentially 
independent of T. The above equation can be used to give a crude estimate of n n 
required for an r process: 
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where nominal values for (vo n j(Z,A)) and A«(Z,A) have been used. 

As a parcel of material undergoes r-processing, its n n and T evolve with time t. The 
temperature T is crucial in determining the path of the r process. If T is always far 
below 10 9 K, the r process is cold and photo-disintegration is negligible. In this case, 
the r-process flow is controlled by the competition between neutron capture and /3 decay, 
and the path is populated dominantly by nuclei with Xp(Z,A) ~ n n (vo n _y(Z,A)) . 

If T significantly exceeds 10 9 K initially, a hot r process occurs. For T ~ 2 x 10 9 K and 
n n 10 20 cm~ 3 , photo-disintegration is sufficiently fast to establish an equilibrium with 
neutron capture. In this (n, y) ^ (y,«) equilibrium, the total abundance of nuclei with a 
specific Z is concentrated in a few so-called waiting-point nuclei, which are favored by 
the relevant n n and T based mainly on their neutron separation energies S n {Z, A) and to a 
lesser extent, on their nuclear partition functions. As with any statistical equilibrium, the 
exact rates do not matter so long as they are fast enough for the equilibrium to hold. So 
the r-process path formed by the waiting-point nuclei is specified by n n , T, and S n (Z,A) 
but is insensitive to (v<7 W) y(Z,A)) and Ay ;M (Z,A) when (n,y) # (y,n) equilibrium holds. 
The /3 decay of the waiting-point nuclei regulates the relative abundances for different 
Z and dominantly controls the progress of the r-process flow. However, when n n and T 
drop below some critical values, (n, y) ^ (y,«) equilibrium breaks down. This allows 
(v(7„ jr (Z,A)) and Ay ; „(Z,A) to play important roles in determining the final abundance 
pattern produced by a hot r process. 

For both cold and hot r processes, the overall abundance pattern produced in an in- 
dividual astrophysical event, e.g., a core-collapse supernova, is a superposition of abun- 
dance patterns produced in different parcels of material that have undergone different 
time evolution of n n and T. Astrophysical models of a hot r process are reviewed in the 
next section, followed by discussion of a cold r process induced by neutrino interaction 
in supernova helium shells. 



ASTROPHYSICAL MODELS OF A HOT r PROCESS 



In the absence of fission, the abundance pattern produced in a parcel of material under- 
going a hot r process can be approximately characterized by the corresponding average 
mass number (A r ) . By mass conservation, 

(A r ) = {A s )+n/s, (3) 

where (A s ) is the average mass number of seed nuclei that capture neutrons and n/s is 
the number ratio of neutrons to seed nuclei at the beginning of the r process. So the 
conditions required for a hot r process to produce nuclei with an average mass number 
(A r ) can be mapped to those giving rise to the corresponding (A s ) and n/s. The viability 
of a model for a hot r process can then be assessed by how seed nuclei are produced and 
how the neutron-to-seed ratio n/s is determined before the actual r process occurs. 

A generic sequence of evolution prior to a hot r process can be illustrated by consider- 
ing a parcel of material with an entropy S > 10 in units of Boltzmann constant per baryon 
that adiabatically expands from T ^> 10 10 K. The material at T J> 10 10 K is essentially 
composed of free nucleons with Y n « 1 — Y e , where Y„ is the number fraction of neutrons 
and Y e is the electron fraction (net number of electrons per nucleon). Given S and Y e , the 
nucleosynthesis during the expansion can be calculated once the time evolution of T is 
specified. A convenient parametric form is T(t) = T(0) exp(— t/Zd yn ), where Td yn is a 
constant dynamic timescale. 

For Y e < 0.5, nearly all protons are assembled into 4 He when T decreases to ~ 
5 x 10 9 K. At lower T, an a process occurs to burn neutrons and 4 He [8]. Here the 
role of S is manifest. Due to the absence of stable or long-lived nuclei at A = 5 and 
8, the nuclear flow must go through the three-body reaction 4 He + 4 He + n — > 9 Be + y 
to produce nuclei beyond 4 He. For S > 100, there are a sufficient number of photons 
at or above 1.573 MeV, the threshold to break up 9 Be, in the Planck distribution for 
5 x 10 9 J> T J> 3 x 10 9 K. So the higher S is, the more inefficiently heavy nuclei are 
produced. In any case, charged-particle reactions effectively cease when T drops to 
~ 3 x 10 9 K due to the Coulomb barrier. By then the inventory of heavy nuclei produced 
typically have A ~ 60-90. These nuclei act as seed nuclei to capture the remaining 
neutrons during the subsequent r process as T continues to decrease. 

Clearly, there are an infinite number of combinations of S, Y e , and Td yn that can give 
the same (A r ) according to Eq. (3). This can be seen from the determination of n/s. 
Approximately speaking, there are three different ways to produce a higher n/s: (1) a 
lower Y e to increase the initial neutron abundance, (2) a higher S to suppress production 
of seed nuclei by keeping the abundance of 9 Be low, and (3) a smaller Td yn also to 
reduce seed production but by shortening the duration of the a process. The r-process 
abundance pattern in the solar system is characterized by two peaks at A ~ 130 and 195, 
respectively. To produce the peak at A ~ 195 requires (A r ) 200. Figure 1 shows three 
contours in the (S, Y^-plane that would all give (A r ) « 200 for the corresponding values 
of Tdyn indicated. These results are taken from [9] and can be used to gauge the viability 
of an astrophysical environment to produce an r-process abundance peak at A ~ 195. 
Similar results can be found in [10, 1 1]. Specific astrophysical models of a hot r process 
are discussed below. 




0.5 - 



I I I I I I I i i u I I I I I I i i I I I I I I I I i i I I I I I I I I i i I I I I I I I I i i I I I I I I I I I I 

100 200 300 400 500 6 



i i i i i i i i i | 1 1 i i i i i i i | i i i i i i i i i | i i i i i i i i i | i i i i i i i i i | i i i i i i i i i 



600 



S 



FIGURE 1. Combinations of S, Y e , and Td yn required for adiabatically expanding material to produce 
nuclei with (A r ) « 200 through a hot r process. See text and [9] for details. 



Core-collapse supernovae are intimately associated with neutrinos. Essentially all 
the gravitational binding energy of a proto-neutron star formed in these events, ~ 
3 x 10 53 erg, is radiated in v e , v e , v M , v M , V T , and v T over ~ 10 s. A fraction of the 
v e and v e can be absorbed by the material immediately above the proto-neutron star 
through the reactions 



During the ~ 10 s of neutrino emission, the proto-neutron star blows a wind as the energy 
deposited by the above reactions continuously lifts shells of material out of its gravita- 
tional potential well. If reaction (5) proceeds more efficiently than reaction (4), the wind 
material is driven neutron rich (Y e < 0.5) [12]. When this material reaches sufficiently 
large radii, these reactions effectively cease due to small neutrino fluxes. Subsequently, 
the material approximately undergoes the generic evolution outlined above for adiabatic 
expansion with fixed S, Y e , and Td yn - 

Earlier studies of nucleosynthesis in neutrino-driven winds based on the supernova 
models of J. R. Wilson showed that a hot r process occurs (e.g., [13, 14]). However, 
later studies [15, 16, 17, 18, 19, 20] could not reproduce the high entropies in these 
models. Typically, 10 £ 100, 0.4 ^Y e £0.5, and 0.01 & T dyn < 0.1 s were obtained. 
For these conditions, the elements Sr, Y, and Zr with A ~ 90 are readily produced, and for 
the most favorable among these conditions, the elements from Pd and Ag with A ~ 110 
up to Te immediately below A ~ 130 can be produced (e.g., [9]). However, all these 
elements are produced in quasi-equilibrium of forward and reverse reactions involving 
photon, neutron, proton, and a particle but not by an r process (e.g., [8]). In fact, heavy 



Neutrino-driven winds from proto-neutron stars 



v e + n ->■ p + e , 
v e + P -> n + e + . 



(4) 
(5) 



r-process elements with A > 130, especially those in the solar r-process abundance peak 
at A ~ 195, cannot be produced under these conditions (see Fig. 1). 

As Y e depends on the competition between reactions (4) and (5), it is extremely 
sensitive to the difference between the emission characteristics of v e and v e [12]. In 
contrast, S and Td yn are not sensitive to this difference as they depend on the net heating 
by both v e and v e . Recent calculations including more neutrino transport processes 
showed that neutrino-driven winds are always proton rich (Y e > 0.5) [19, 21]. This 
gives rise to a new process of nucleosynthesis, the vp process [22, 23], which produces 
elements beyond the Fe group mostly via proton capture. Without neutrinos, those 
"waiting-point" nuclei in (p, y) ^ (y,p) equilibrium must /3 decay before the nuclear 
flow can proceed further (cf. a hot r process) and a waiting-point nucleus with a very 
slow /3 -decay rate effectively terminates the flow. The flow proceeds much further in the 
presence of an intense v e flux because the neutron produced by reaction (5) can enable 
an (n,p) reaction that has the same effect as but is much faster than /3 decay. Therefore, 
the vp process can produce elements far beyond the Fe group. It was shown that the 
elements from Sr, Y, and Zr to Pd and Ag with A ~ 90-1 10 can be produced in both 
proton-rich and neutron-rich winds [24]. So while neutrino-driven winds are unlikely 
to produce heavy r-process elements with A > 130, they are an important source of 
the above elements. Further, they served as a detailed case study of novel processes of 
nucleosynthesis that were not envisioned in the classical papers of Burbidge et al. [1] 
and Cameron [2] on the origin of the elements. 



Shocked surface layers of O-Ne-Mg cores 

Stars with masses M J> 1 1M© develop Fe cores and those with 8 %,M < 1 1M© develop 
O-Ne-Mg cores before they undergo core collapse (e.g., [25]). These two groups of stars 
differ greatly in structure as the density above an O-Ne-Mg core falls off much more 
drastically than that above an Fe core [26]. So far O-Ne-Mg core-collapse supernovae 
are the only case for which neutrino heating [mostly via reactions (4) and (5) as in the 
case of neutrino-driven winds] has been demonstrated to result in a successful explosion 
(e.g. [27]). 

The surface layers of an O-Ne-Mg core are slightly neutron rich due to processing 
by pre-supernova nuclear burning. The evolution of such a layer subsequent to being 
struck by the supernova shock again approximately follows adiabatic expansion with 
fixed S, Y e , and Td yn . Due to the steep density gradient above the O-Ne-Mg core, the 
shocked material may expand rapidly with Td yn ~ 10~ 3 s. This along with S ~ 100 and 
Y e < 0.5 (see Fig. 1) would enable a hot r process to produce nuclei with (A r ) ~ 200 
[28]. However, current models of O-Ne-Mg core collapse based on Nomoto's models 
of the associated supernova progenitors [26] only give a weak explosion, for which the 
shocked surface layers of the core expand with Td yn 3> 10~ 3 s and S <C 100 [29]. As these 
layers are not extremely neutron rich, no hot r process can occur under the conditions 
obtained in current models. It remains to be seen if more favorable conditions can be 
obtained for new models of stars that develop O-Ne-Mg cores. 



Winds from accretion disks of black holes 



Figure 1 shows that for T dyn /S 3 ~2x 1(T 7 s with 25 < S < 100 and 3.9 x 10~ 3 < 
T dyn ^ 0. 195 s, nuclei with (A r ) 200 can be produced by a hot r process for Y e ~ 0.2. 
The above ranges of S and Td yn occur in neutrino-driven winds from proto-neutron stars. 
However, Y e ~ 0.2 cannot be obtained in these winds without invoking new physics 
such as sterile neutrinos [30]. Therefore, it is important to explore alternative candidate 
environments for a hot r process where such low Y e values may be obtained more 
naturally along with similar S and Td yn to those in neutrino-driven winds. 

Some Fe-core collapse supernovae are expected to produce black holes rather than 
neutron stars. An accretion disk typically forms around the black hole. Such disks can 
also form in mergers of two neutron stars or of a neutron star with a black hole. The 
disk material close to the black hole is hot and dense and emits neutrinos similarly to 
the surface layers of a proto-neutron star. In addition, this material may be very neutron 
rich as a result of the reaction e~ + p — > n + v e . A fraction of this material is ejected 
in winds due to heating by dynamic processes in the disk or by neutrino reactions (4) 
and (5). Parametric models suggest that adequate conditions, in particular low Y e values, 
may be obtained in these winds for a hot r process to occur [31, 32, 33, 34]. It remains 
to be seen if such conditions can be derived in self-consistent models of the accretion 
disk and the associated winds. It is especially important that detailed neutrino transport 
in the disk be carried out similarly to that in modern core-collapse supernova models. 



Neutron star mergers 

Material of low Y e stocked up inside a neutron star can be released by tidal force 
during mergers of two neutron stars or of a neutron star with a black hole. Once released 
from the strong gravitational hold of its parent neutron star, such material undergoes 
decompression, which is accompanied by heating due to /3 decay. The possibility of a 
hot r process during the expansion of decompressed neutron-star matter was discussed in 
the pioneering work of Lattimer et al. [35] and more recently in [36, 37]. An interesting 
feature of this r-process model is that the nuclear flow typically reaches fissile nuclei 
and becomes cyclic as the fission fragments capture neutrons to replenish the flow. This 
fission cycling is expected to produce a robust abundance pattern at A > 130, which 
appears to be in accord with observations of metal-poor stars in the Galactic halo but 
with some notable exceptions (e.g., [38, 39]). 

A potential problem with r-process models based on neutron star mergers is that 
these events occur ~ 10 3 times less frequently than core-collapse supernovae in the 
Galaxy. Were neutron star mergers a major source for heavy r-process elements with 
A > 130, stars with the earliest such enrichments would have been enriched with high Fe 
abundances by the numerous core-collapse supernovae that had already occurred prior 
to the contributing r-process event [40, 41]. This is in contradiction to observations of 
a number of metal-poor stars that are enriched in heavy r-process elements but have 
[Fe/H] = log(Fe/H) — log(Fe/H) < —3 (e.g., [38]). It remains to be seen if more 
detailed models of Galactic chemical evolution may resolve the above difficulty. 



NEUTRINO-INDUCED NUCLEOSYNTHESIS IN SUPERNOVA 
HELIUM SHELLS: A COLD r PROCESS 



The notion of a unique type of sources that produce heavy r-process elements with 
A > 130 over the entire history of the Galaxy has some appeal but is far from required 
by observations (e.g., [39, 42]). As discussed above, winds from black hole accretion 
disks and decompressed ejecta from neutron star mergers could both be production 
sites for such elements. The difficulty of neutron star mergers to account for observed 
enrichments of heavy r-process elements at [Fe/H] < —3 leads to the question whether 
other sources could have operated at such low metallicities. In fact, such a source was 
proposed by Epstein, Colgate, and Haxton (ECH) [43]. In their model, NC neutrino 
reactions on 4 He nuclei produce neutrons for an r process via 4 He(v, vn) 3 He(n, p) 3 H or 
4 He(v, vp) 3 U followed by 3 H( 3 H,2n) 4 He. This occurs in helium shells of metal-poor 
core-collapse supernovae, where neutrons are left alone by the predominant 4 He nuclei 
but captured by the Fe nuclei acquired at the birth of the supernova progenitors. The 
helium shells have T ~ 10 8 K, and therefore, the ECH mechanism is a cold r process. 
The neutron density for this mechanism depends on the competition between production 
by neutrino reactions and capture by Fe and other nuclei. Consequently, n n J> 10 18 cm" 3 
required for an r process [see Eq. (2)] can only be obtained at low metallicities. The ECH 
mechanism was critiqued in [44] as viable only for low-mass metal-poor supernovae, 
where helium shells are located at relatively small radii to ensure sufficiently large 
neutrino fluxes. Follow-up studies were carried out in e.g., [45]. 

Recently, it was shown that a modified ECH mechanism could operate in core- 
collapse supernovae with [Fe/H] <, —3 over a range of progenitor masses [46]. For a 
\5Mq progenitor with [Fe/H] ~ —4.5 [25], a cold r process occurs in the outer helium 
shells at radii of ~ 10 10 cm, where the abundances of 12 C and 16 nuclei produced by 
pre-supernova helium burning are too low to cause significant neutron capture. The CC 
reaction 4 He(v e , e + n) 3 H turns out to be the dominant source of neutrons because the 
temperature is too low for the reaction 3 H( 3 H,2«) 4 He to be effective [46]. Further, the 
CC reaction rate is greatly enhanced when the energy spectra of and v T are much 
harder than that of V e at emission and an inverted neutrino mass hierarchy enables 
v e Vji.T transformation below the outer helium shells [46]. More systematic studies of 
the modified ECH mechanism are needed to assess its contribution to Galactic chemical 
evolution. It is conceivable that this mechanism operated in the early Galaxy until 
[Fe/H] ~ —3 and neutron star mergers provided the subsequent r-process enrichments. 
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